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DERIVATION OF FORMULAS FOR EVALUATING THE STANDARD 


ERRORS IN B, C, AND N,_) WHICH APPEAR IN THE 
r 


EQUATION: Z_ = (Z./P,)£ Noug P. 


by 


1/ 


B. J. Dalton and Robert E, Bares 


ABSTRACT 

The method of treating a set of isothermally measured pres- 
sures Po? Pi Po> we tS Pa? for a Burnett experiment, consists of 
expressing the compressibility factor isotherm of the gas in terms 
of a function of either P or p and evaluating the volume ratio, 
Np=0° and the constants in the function by least squares solution. 
This report gives the method of least squares as applied to the 
fundamental equation for the Burnett experiment, taking into con- 


sideration the change in N with pressure. 


P=0 
Formulas are given for evaluating the standard deviation of 


a single Pe the standard deviation of each of the constants eval- 


uated, and the standard deviation of the compressibility factor. 


1/ Research Chemist, Helium Research Center, Bureau of Mines, 
Amarillo, Texas. 
2/ Supervisory Research Chemist, Project Leader, Thermodynamics, 


Helium Research Center, Bureau of Mines, Amarillo, Texas. 
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INTRODUCTION 
The method of treating Burnett type data consists of expressing 
the compressibility factor isotherm of the gas in terms of a func- 
tion of P or p and evaluating the best values for the constants in 
the function by least squares solution. 


In a previous report Gee we have given formulas for eval- 


3/ Underlined numbers in parentheses refer to items in the list of 


references at the end of this report. 


uating the best values for the constants appearing in the fundamental 
equation for the Burnett experiment: Z = (Z./P,)N as assuming N, 
the cell constant, to be independent of the pressure. 

Canfield (5, 7) has pointed out that the volume ratio varies 
with pressure due to distortion of the bombs and a shift in the null 
point of the differential pressure indicator. This change in the 


cell constant, for a given expansion, was expressed by the equation 


IL+aqaP. 
Nee = Nee yee (1) 
r P=0\ 1 3°68 Peel 
where 
N5-0 = volume ratio at zero pressure 
N.. = volume ratio for the ae expansion 
th : 
P = pressure after the r— expansion 
I = pressure before the ae expansion 
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@P = change in the total volume, 
WV, + V2) Pep divided by the 
total volume at zero pressure 

BP oy = change in vy at P=P ol divided 
vy at P=0 


Therefore, representing the change in the volume ratio by equation 


(1), the fundamental equation for the Burnett experiment assumes 


the form 
Z = (Z /P .)N > N, ; NP. (2) 
Since 
. " (ee cat + ge sou) 
1 ag O\L + BP. 
: € + 2) 
Ps P=0 \l + BP, 
1 + oP 
N = N ( + ) 
r P Lor 6P 
r 
equation (2) is expressible as 
: r 
Z = (Z /P DE N50 B (3) 
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Associated with this function, equation (3), whose constants 
have been evaluated by least squares, the following errors are of 


interest: the standard deviation of a single P the stan- 


r (obs)? 
dard deviation of each of the constants; and the standard deviation 
of compressibility factors. 

The method outlined in reference 1 served as the basis for 
evaluating the least squares solution for the constants appearing 
in equation (3). The methods outlined in reference 2 were used 
for evaluating the above-mentioned errors. 


METHOD OF OBTAINING THE LEAST SQUARES SOLUTION FOR THE 
CONSTANTS APPEARING IN THE EQUATION: Z = (Z /P Et N5=0 ie 


The fundamental equation for the Burnett experiment, assuming 
the variation of the volume ratio with pressure to be expressible 


by equation (1), is of the form 


r 
Peo aiat Oe ee ay alana (3) 


where £ is defined by equation (4). 

Now there are two series expansion which can be employed for 
representing the compressibility factor isotherm: the Leiden 
expansion in powers of p and the Berlin expansion in powers of P. 


The Berlin expansion 
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was chosen because all of the parameters for which we seek a least 
Squares solution can be expressed in terms of the original obser- 
vations. It was assumed that the series expansion in powers of P 
could be truncated after the third virial coefficient. 

Suppose we let the functional relationship between the variables, 


P. and r, involving the three parameters, B, C, be 


Np=0? 


FE. =, EGu,CPs> N 


r P=(0" B; C) bs 0 (7) 


Now because of random errors in the observed pressures, when yee) 
is substituted in the above expression, F will not be exactly zero. 


Let F. be the value of F when the observed values of r and P. are 


substituted in equation (7). Thus, 


EF. wen ir, B, C) (8) 


r Fe (obs) ? Np=0? 
Now we assume that r, the expansion number, is accurately known. 
Therefore, equation (7) may be solved for Shc) so that equa- 


tion (7) is exactly satisfied. Thus, 


F = F(r, B, C) = 0 (9) 


Beale)’ Np=0? 


and P is the solution of equation (9). 


r(calc) 
Now AP > the residual of Pe? is the difference between the 
observed and calculated values. This is not the true random error 


in our observed P. because we do not know the true value of hs, 


However, we can maximize the probability that the AP's are equal 


ord 1) G a4 i. 7 
tence! » dose ww Hokie wi laa sda to. ie et nee wots ase 


~teedu fachaktq sit to mri29 " Sentence od aa oka utoe isi 7 


7 lo avswog nt colensyes eatiee off tad9 bemswes caw 3D ,en0k se: er 
- é rn as ry 
.Insiottiacea fatuty betdo silt asste baceomisd sd binds 


,eaidsatisy 249 asswied ¢cidenoisales Jenoitjonn) af3 ial aw Sesgqua 


o9 .9 .a4, ¥_. ered sq Stl ja ativioun eet 
9 Fl egg maszeq somid sia ativiounti .x bus 3 
{<) f} = f " j , g 5) = q 7 
4 (= ; —_ 4 
4 ™ 3 owtiw .2Siudeang bavieedo sfi at @26555 mobaet to sdunced wok 
4 
’ a E — 
oz9S yvijosxs sd jon iltw 71 _,u9lméeages svoda eds at beausizgedue ai 
2 q cs 30 aSsHley 6 ‘Seno 3n3 nearit i 3o oul si} od i 4a! 
. aud (jy) Geis Haps af bagusisedue e 
7 
: ) ) : 4 - ~~ > @ Ba Und 
. ‘o=q" *(edo)2 +* we 
4 viestsmwoos ef ,Istdawi doltensexre St ,1 Jsadd ameesa ow woul 
r, a 
: . . . , | 
indd of ‘ Y tot bavioe ad vam {(V) aolisaups ,syolseta!l ' 
(ofao)a . mo 
, ,auaqT .botleitae via xsxs st (4) nors 
e) O a= 0 fs » — .%)7' = . 
™~ ©" “tT ‘Colés)y * => 
1 fg 4 Les - 
.(@) wofiawps io gotsyviog sid at q bos . 


(o489)% : 7 


7 = 
ij geswisd soas7S)ikb ond af , 7 to Inubiesy ols , a Wo 
7 Oe . 
ei 


a 
mak 


7017x989 mobras oia9 oft Bon ob aia -2oulav saan a rime _ vary 


-,i Jo sufsv sute od? hers ah wes - 
: ors Fata! 
wps ois a’ x a te dnd is my ote : bn 


ad 


. ple 


a 


to the true random errors, and this is just what the principle of 
least squares does. The principle of least squares says that we 

maximize the probability that the AP'S represent the true random 
errors by minimizing the sum of the squares of the weighted resi- 


duals. Thus, we should minimize the function 


ie 
R= ) W (ap)? (10) 
af r (obs) 


and evaluate N Beegand C so that 
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then they will have the same weight and W a) boo phe Pos 
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do not all have the same precision index, then 
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In a particular problem, it may be necessary to assume 


Wp = 1 in the beginning. However, if this is done the 
r (obs) 


residuals, Y. = LE (aba) = Pirate: should be examined to see 
if there is any statistical evidence for the residuals squared 


being a function of P Any assumption as to the variance 


r(obs)~ 


being a function of P can always be checked by examining the 


r (obs) 
residuals. In any event, W is not a function of the con- 
r (obs) 
stants to be evaluated. 
In order to evaluate Np-o° B, C, we need to linearize Y, 
with respect to the undetermined constants. A truncated Taylor's 
series expansion was used to do this. 


In a previous report (1), we show that the linearized normal 


equations are expressible as 


a, AB + b, Ac + c ANZ _9 = mm (14) 
a, AB + b, Ac + c, AN5_o = ml, (15) 
ajhB+ bsACe+ e,AN,_,) = mm, (16) 


Equations (14), (15), and (16) result from expanding Y> (o¥, /OB) , 
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where B, C, N5-9 are the undetermined constants and B ,C , N5-0 


are approximate values for these quantities. 
The a's, b's, c's, and m's appearing in the normal equations 
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Now in evaluating the best values for B, C, and N =o we solve 


i 
the linearized normal equations by an iterative procedure. Even 
though it is necessary to solve equations (14), (15), and (16) by 
a series of approximations, it is important to realize that the 
best values for the constants are determined such that the normal 


equations are exactly satisfied. 


The solutions to equations (14), (15), and (16) are (1): 


D AB = Dim, + Dom, ay Dm. (39) 

D Ac = Dm, + Dem, + Dem, (40) 

D AND. = Dom, + Dem, + Dom. (41) 
where 

D, = boc, - b,c, (42) 

D, = D, = bac) - bic, (43) 

MES cate D1 2 Oot or 

D,. = ajCg - ancy (45) 

De = De = a5) - aslo (46) 

Dy = a,b, - ayb, (47) 

D = Dia, + Daa, + Daa (48) 
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EXPRESSIONS FOR DETERMINING VARIANCES 
AND COVARIANCES OF THE CONSTANTS EVALUATED 


We now proceed to evaluate the variances and covariances of 
the constants evaluated. To do this, we apply the law for the 
"Propagation of Errors" (3, 8). This law states that if we have 
a function or quantity, say Q, that is a function of the indepen- 
dently observed quantities Yy> ¥5s oes then the variance of the 


quantity Q is given as 
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Ze : Z : : 
where Sq is the variance of Q and §S is the variance of 
Yi (obs) 


Yi(obs)’ Extracting the square root of the variance, we obtain 


a value on the same scale as the original measurements. This 


value, S., is called the standard error or the standard deviation 
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Since we have assumed that r is accurately known, then the expres- 
sion for the variance in B is given by the equation 
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and there will be an equation similar to equation (50) for deter- 
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In order to evaluate equation (50), we must evaluate 


(OB/dP ) for each P )? multiply this quantity by Sp : 


ECobs r (obs) 


square the product, and then sum the product over all of the 


rt (obs 


observed PS 


In a previous report (2), we have outlined the details for 
evaluating the variances and all of the covariances of the con- 
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We have indicated that the solutions to equations (14), (15), 
and (16) are solved by an iterative procedure. When we have the 
correct solutions, our linearized normal equations will be satis- 
fied exactly. Therefore, once the best values for the constants 
have been determined, the remaining questions to be answered are: 
(1) what is the variance of the calculated P's and any other cal- 
culated P that reduces F to zero?; and (2) what is the variance 


of the compressibility factor? 
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where 
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EXPRESSION FOR EVALUATING THE VARIANCE OF THE COMPRESSIBILITY FACTOR 
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Now in order to evaluate equation (68), we must evaluate 
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ee eaic). robs)! can be determined from equation (59), 

where the terms involving derivatives of F are to be evaluated 
; £ e e e 

for Ficale)’ rp fp, and [OF/OP ) sic)! is defined by equation 


(67). Multiplying equation (59) by (aB/dP ie and 
r (obs) E (obs) 


summing the product over all of the observed PSs we get 


a 
ce 
B\OB 
| ee re0" (calc) 
2a oF? 
+ S3¢ (SE). Nt see 
Poo PeO7 (calc) 
2 ( OF ) 
E ,aP i: BENee aU pap rent cer 
i (calc)\/ 8B \ 2 oe? J pB,C,P (c 
cae SE = § = = (72) 
Z, \OP oP P B,P OF 
r (obs) r(obs) r (obs) (calc) a 
— pErenic ro Sb Gl 
Pp???" p= 
; 2 
If we multiply equation (59) by (8C/OP (obs) Sp and then 


| r (obs) 
sum the product over all of the observed PS we get 


BC \oB 
Postel pef* (cate) 
2@ 
Gece resp N P 
P?~? Pa? (calc) 
r he ce 
cate) P=0 ‘“"p=0’ r_,B,C,P 
i peetaba) fobs) aes) Oe (Cealcy =) 
r= (calc) ¥p»B,C,Np_p 


beseulave ad 03 916 4 20 | 
notsaups yd baniteb at eos tis eos 
bre rCy, foystONAS) ed (RA) nokteups g 


cada? 
jog 30 e' 4 bevisedo of? to Tin xsvo tae 


a 


(shes) "00g" 


~ 
a a ca 


Sp 
(eda) z { edo) nie 


angle 9s Bn gt (oles) 


, 
ned? bem ( T6\96) yd (RC) wolseups vigislum ow 22 
a woud (edo) 


jsg sw ert beyteedo efi to Ils xevo Jouborg sii ma | 


(36 ¢ . 
@ | 3 
si A) oa | . | 


(ofma) 2 Qaq@s rg 


36 


Therefore, the use of equations (72) and (73) in equation 


(69) will enable the variance of a calculated Zp» sf rated 
Brealey 
be determined. 
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